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Abstract 

A decaying particle that dominates the energy density of the universe and then 
decays with a lifetime ~ 10"^ sec into electromagnetically-interacting particles can 
heat the background electron neutrinos more than the fi and r neutrinos, resulting 
in a decrease in the primordial ^He abundance. Previous calculations of this effect 
assumed sudden decoupling of the neutrinos from the electromagnetic plasma. Here 
we repeat this calculation using an improved (although still approximate) treatment 
of the energy flow from the blackbody radiation into the various neutrino components. 
We find that with our improved calculation, the reduction in the primordial ^He is 
not as large as had been previously thought, but it occurs over a larger range in 
decaying particle lifetime r. The reduction in ^He can be as large as AY = —0.01 for 
r = 0.1 -0.7 sec. 
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I. INTRODUCTION 

The effects of decaying particles on primordial nucleosynthesis have been inves- 
tigated by many authors (see, for example, [1] - [2] and references therein). One 
interesting effect, first noted by Scherrer & Turner [3] can occur if the particle de- 
cays into products which interact electromagnetically and go into thermal equilibrium 
with the cosmic blackbody background. If such decays occur well after the neutrinos 
have dropped out of thermal equilibrium (at a temperature of a few MeV) then the 
neutrinos do not share in the heating of the blackbody background. For hfetimes 
much shorter than the neutrino decoupling temperature, the neutrinos share fully in 
the heating of the blackbody background. In the intermediate regime (lifetimes on 
the order of 0.1 sec), an interesting effect can occur. Because the electron neutrino 
couples more strongly to the e~e'^ pairs than do the /i and r neutrinos, it decouples 
at a slightly lower temperature. Thus, for lifetimes in this intermediate range, it is 
possible for the decaying particle to heat the electron neutrinos more than the /i and 
T neutrinos. This results in a net decrease in the expansion rate relative to the weak 
{n <-> p) rates, giving a reduction in the ^He abundance [3] . This effect can occur, for 
example, in gravitino decay [4]. 

This possibility of reducing the primordial ^He production is quite interesting, be- 
cause recent calculations suggest that the standard model of big bang nucleosynthesis 
can be made consistent with current observations only if the "true" primordial ^He 
abundance is larger than currently believed, or if primordial nucleosynthesis is modi- 
fied to give a smaller production of ^He [5]. However, the treatment of this effect in 
references [3]- [4] was quite crude: the neutrinos were assumed to decouple suddenly 
when their interaction rates dropped below the expansion rate; entropy dumped into 
the system before decoupling was shared equally by the neutrinos, photons, and e'~e'^ 
pairs, while entropy released after decoupling went entirely into the photons and e^e"*" 
pairs. Hence, we are motivated to do a more careful investigation, involving a more 
accurate treatment of the transfer of energy from the decaying particle to the neu- 
trinos. We use the linearized energy-transfer equation given by Rana and Mitra [6], 
based on earher work by Herrera and Hacyan [7] . We do not calculate the full distor- 
tion in the spectrum, but instead treat each species of neutrino as a black body with 
a single temperature. Although this is not as sophisticated as the treatment given 
by Dodelson et al. [2], it represents a significant improvement over other previous 
studies of this problem. (Dodelson et al. examined only decaying neutrinos, which 
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do not have a large enough energy density to make this an interesting effect) . 

In the next section, we present our calculation of the evolution of the neutrino 
temperature in the presence of a decaying particle with a short (10~^— 2 sec) hfetime. 
In Sec. Ill, we calculate the effect of the differential neutrino heating on the primordial 
^He abundance, and briefly summarize our conclusions. We find that the reduction 
in ^He production is smaller than had been previously thought, but it is still large 
enough to be interesting. 

II. EVOLUTION OF THE NEUTRINO TEMPERATURE 

The calculation begins by adding to the standard cosmological model an X particle 
which decays exclusively into particles that thermalize rapidly compared with the 
expansion timescale. This particle adds a contribution to the total energy density 
driving the expansion according to the Friedmann equation 

H^^ = {lnGpy/'. (1) 

where R is the expansion scale factor. Here the total density p = pstd + Px with psta 
representing the energy density in the standard model. The X density px evolves 
according to the equation 

Px = -rpx - 3i/px, (2) 

where V is the decay rate, T = l/r, and r is the X particle lifetime. This equation 
has the solution 

PX = Pxoi^r'e-^'^. (3) 

We consider particle hfetimes in the range 10~^sec < r < 2 sec. Following reference 
[3], we quantify the the particle number density by the parameter r, which gives the 
number density of the X particles relative to photons before e~e^ annihilation. The 
only quantity which determines the energy density of the X particles as a function 
of temperature is rmx, where mx is the particle mass. In all of our calculations, 
we take rmx — 10^ MeV. For this choice of rmx and particle hfetimes above, the 
density of the X particles dominates the expansion prior to decay. More importantly, 
the entropy produced by the decay dominates the previously-existing radiation during 
nucleosynthesis. When these two conditions are satisfied, the evolution of the element 
abundances is independent of rmx, and the final results will be essentially a function 
only of T [3]. For the range of hfetimes considered here, these conditions will be 
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satisfied for rmx > 10^ MeV, so tliis is tlie range over wliicli our results will be 
applicable (See ref. [3]). For comparison, a massive tau neutrino has a maximum 
rnix of about 1 MeV (for m,y ~ 3 — 5 MeV), so we are examining a very different 
region of parameter space from that discussed in ref. [2]. 

The decay of the X particle heats the electromagnetic plasma. Let and be 
the energy density and pressure of the photons and all relativistic particles in thermal 
equilibrium with the photons. Then the evolution of is given by 



In references [3]- [4], the authors used a sudden decoupling approximation to deal 
with neutrino heating. In this approximation the neutrinos are assumed to have 
blackbody spectra with temperatures T^,-, i = e,/i, r. The temperature of an indi- 
vidual neutrino species is set equal to the photon temperature until it decouples and 
thereafter the neutrino temperature decreases as 1/R. The decoupling condition used 
in references [3]-[4] is T^JH < 1 where Vy- = A^T^ {i = e,/i, r). For consistency 
with our derivations below, we re-express Ti,. in the form 



where is the Fermi coupling constant, and a^- is a dimensionless constant of 
order unity. The values for A^- used in references [3]- [4] correspond to a^^ = 0.3 and 
a^^ ^ = 0.07. This is obviously a very crude approximation and is considered here 
only for comparison with our new calculations. 

In this work we follow the treatment of Rana and Mitra [6] and Herrera and 
Hacyan [7], who treated the neutrinos as a perfect blackbody and calculated the rate 
of energy transfer from the electron-positron pairs to the neutrinos through scattering 
and annihilation. Let u^. be the total rate of energy density transfer from the e~e"'" 
pairs to a single neutrino species z/j [i = e, p, r) via annihilations and scatterings. 
Then equation (^) is modified to 




(4) 



(5) 



p^ = Tpx 



Y: u,^-3H{p + p) 
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(6) 



while the neutrino density evolves as 



= - AHp^., 



(7) 
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In the limit of small temperature difference between the electrons and neutrinos, Uy^ 
is given by [6] 

= 2Je-^,,(T^-T,J (8) 

where le-^ui is a function of which gives the rate of energy density transfer per 
unit temperature difference from electrons into a particular neutrino species i>i, and 
the factor of 2 comes from including both electrons and positrons. In the limit of 
large temperatures (T^ ^ me), / has the form 

h-^Ui = Pe-^UiG'pT^, (9) 

where jSe-^ui is a dimensionless constant of order unity. [Equation comes from the 
fact that le-^ui is approximately equal to the interaction rate (F = n{av) ~ G'^pT^) 
times the heat capacity (~ T^).] Three distinct processes contribute to the net value 
of /e-^j,-: e~e^ annihilation into ViVi pairs (along with the inverse annihilation), e~ 
scattering off of z/j, and e~ scattering off of z/j. The values for (3f,-^y- can be derived 
from the high-temperature limit of the results of Rana and Mitra [6] and are listed in 
Table 1, where we have included both the total Pe-^ui used in equation (H), as well 
as the contribution to f3e-^ui from the individual reactions. It is obvious from the 
numbers in Table 1 that e'e'^ annihilation is the dominant energy transfer mechanism. 

In calculating the relative temperatures of Ue, z/^, and Ur we must also consider 
the rate of energy transfer from the Ue into the other two neutrino types, a process 
neglected in reference [6] . Again, we can model these rates in terms of an energy loss 
term w^-, i = /i, r, where 

w,^ = 2h^^,XTu.-nj (10) 
and lue^Vi is a function of T^^ which can be written in the form 

^Ue^Ui = Pue-^UiGpT^^ (11) 

Again, there are three processes which contribute to the net energy transfer from 
z/e into uf. z/gZ/e ^ z^jZ^j, z/gZ/j — * z/eZ/j, and UgUi —y z/gZ/j. Using the standard matrix 
elements (see, e.g., reference [8]), we have calculated the differential cross sections for 
these three processes. These cross sections can then be substituted into the integral 
expressions for / corresponding to those given in reference [6] to derive the values for 
Pue^Ui in equation (|TT]). These values are given in Table 1. Again, we see that z/gt'e 
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Table 1: The coefficients for the energy transfer rates in equations (Q) and ([TT]), where z/j 
denotes a or r neutrino. We have included the breakdown of Pe-^u^ iiito specific reactions 
as well as the total Pe-^ui used in equations and (0). 
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7/ 4. /O 7/ 


n 14 

yj . ±t: 


p~ — > p~z7« 


0.04 


e~e^ <— > Uei^e 


0.99 


e~ z/e 


1.17 


e~z/j — > e~z/j 


0.02 


e~z/j ^ e^z/j 


0.02 


e~e~^ ViVi 


0.21 


e~ — > z/j 


0.25 


Z/gZ/i ^ Z/eZ/j 


0.07 


Z/gZ/i VeVi 


0.01 


V^V^ ^ ViVi 


0.42 




0.50 



annihilation is the dominant energy transfer mechanism. With this additional energy 
loss term, equation (|^) becomes: 

Pu, = u^e- - ^Hp^^, (12) 

i = U,T 

Pu, = u^, + - AHp^^ (i = p^t). (13) 

In deriving these equations, we have made two major approximations: first, we 
have neglected the distortion in the energy density spectrum of the neutrinos, treating 
the spectrum as a blackbody and considering only the change in the total energy 
density. This is relatively unimportant for the /i and r neutrinos, since their only 
effect on nucleosynthesis is via their contribution to the total energy density, but it 
could be important for the electron neutrinos, since the weak rates are highly energy 
dependent. Second, equations (||) and (p!0[) are strictly valid only in the limit of 
(Ty — T^.)/Ty <^ 1 and (T^^ — T^J/T^. <^ 1, respectively; we have approximated 
the energy transfer rate by extrapolating these equations into a regime where the 
temperature differences are large. For instance, with r = 1 sec, and for expansion 
times t < 1 sec (the relevant time range for the freeze-out of the n ^ p reactions), 
we find that (T^ - T.^.J/T^ < 0.5, (T, - T^J/T, < 0.35 and (T,, - T,^J/T,^ < 0.25. 
The temperature differences are smaller for shorter lifetimes. For r < 0.1 sec we have 
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m - < 0.2, {T^ - T.J/T; < 0.1 and (T,^ - T.^.J/T,^ < 0.1. So we expect 

our approximation to be quite accurate for r < 0.1 sec, but less accurate for larger 
lifetimes. 

We have integrated our temperature evolution equations for a decaying particle 
with rmx — 10^ MeV. Figure 1 shows the evolution of the neutrino temperatures as a 
function of time for the case of r = 0.1 sec, including the effect of e~e+ annihilation. 
The temperatures are displayed relative to the photon temperature in terms of the 
fraction T^,./T^. We see that for t < t, both T^^^jT^ and Tj,^.^/Try decrease as the 
particle decays, but the electron neutrinos are held closer to the photon temperature 
than the ji and r neutrinos, as in reference [3]. The decrease in T^JT^ which occurs 
at i > 10 sec is the standard effect from the annihilation of the e~e+ pairs and has 
nothing to do with the particle decay. An interesting "rebound" effect occurs between 
10~^ and 1 sec. The neutrino temperatures actually increase briefly relative the photon 
temperature before the e''e^ entropy dump sets in. We believe that this is a real effect, 
not an artifact of the approximations we have used. It arises when the rate of energy 
transfer to the photons (from the decaying particle) is greater than the rate of energy 
transfer from the photons to the neutrinos. Consider, for example, the (unphysical) 
limit where the energy of the decaying particles is transferred instantaneously into 
the photons. After a sharp increase in the photon temperature, the ratio T^/T^ would 
increase as photon energy is transferred to the neutrinos. 

Figure 2 shows the final temperature ratios, multiplied by (11/4)^/^ to factor out 
the effect of e'e^ annihilation, as a function of X particle lifetime. Also shown for 
comparison are the results of the sudden decoupling approximation used in reference 
[3]. This plot demonstrates significant differences between these two treatments. The 
sudden decoupling approximation produces a large difference between the electron and 
other neutrino temperatures in the region r = 0.1 sec. In our more exact treatment, 
the differential heating effect is reduced, but it persists over a larger range in lifetime, 
for T as large as 1 sec. This is what one might have expected, since the more exact 
treatment allows for energy transfer from photons to neutrinos at late times, when the 
sudden decoupling approximation assumes that the neutrinos are already decoupled. 

III. EFFECTS ON ^HE PRODUCTION AND CONCLUSIONS 

We have applied these equations to the primordial nucleosynthesis computer code 
of Wagoner [9], as modified by Kawano [10]. For a decaying particle with energy 
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density rnix = 10^ MeV and lifetime r, we have calculated the change in the pri- 
mordial helium mass fraction Ay as a function of r for three values of the baryon 
to photon ratio: t] = 10~^°, 10~^'^, and 10^^. These results are displayed in Fig. 3. 
We find that AY is nearly independent of r] and can be as large as —0.012. We see 
that Ay < for r < 1.5 sec; for larger lifetimes the decaying X particle, rather than 
its thermalized decay products, dominates the expansion during the freeze-out of the 
n ^ p reactions, resulting in a larger expansion rate and a net increase in '^He. 

For comparison, in Figure 4 we include the Ay value for rj = 10~^'^ using the 
sudden decoupling approximation from reference [3]. We find that our more exact 
treatment produces a smaller decrease in the primordial helium, but that the effect 
occurs over a larger range in particle lifetime; we find that Ay < —0.005 for r = 
0.05—1.2 sec. This is not surprising, since our more exact treatment produces a smaller 
differential heating effect spread out over a larger range in r. Part of the difference, 
however, is that the actual total interaction rates used to compute decoupling in 
references [3] -[4] are different from those used here. So in Figure 4 we also give the 
helium abundances corresponding to sudden decoupling when we change the total 
interaction rates in equation to give the same freeze-out temperatures as those 
derived in ref. [6] {a,^^ = 1.5, au^r — 0.33). The curve is shifted over to correspond 
more closely to our more exact treatment, but it retains a much sharper and deeper 
minimum. 

Of the Ut- decay modes considered by Dodelson et al. [2], our present work most 
closely resembles Ur decay into sterile plus electromagnetic decay products. No re- 
duction in ^He is evident for this case in reference [2], but this is to be expected, 
because the Ut energy density is too low to lead to significant differential neutrino 
heating. Dodelson et al. [2] do see a significant reduction in ^He for decay modes 
which include a u^. in the final state, but this reduction occurs for entirely different 
reasons: the additional decay products increase the n ^ p weak rates, keeping 
them in thermal equilibrium longer. 

Our results indicate that the decrease in ^He production due to differential neutrino 
heating from a decaying particle, first discussed in reference [3], holds up under a 
more exact treatment. The effect is not as large as had been previously estimated, 
but it occurs over a longer range in particle lifetime. The change in Y can be as 
large as —0.01 for r in the range 0.1 — 0.7 sec, and as large as —0.005 for r = 
0.05 — 1.2 sec. A decrease on the order of 0.01 can resolve current problems with 



8 



standard primordial nucleosynthesis [5] , although it is perhaps not the most plausible 
mechanism for resolving these problems. A more detailed treatment of the spectral 
distortion using the full machinery of reference [2] seems justified, although we do not 
expect our results to be significantly altered in such a treatment. Even the treatment 
in reference [2] would need to be modified to treat this problem correctly, since that 
treatment also assumes small temperature differences. 
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FIGURE CAPTIONS 



FIG. 1. The time evolution of Ti,jT^ for electron neutrinos (upper curve) and ji and r 
neutrinos (lower curve) with a decaying particle with lifetime r = 0.1 sec and energy 
density parameter rmx — 10^ MeV. 

FIG. 2. The final ratios {11/ Af'^{T^jT^) function of particle lifetime r for 
a decaying particle with a density parameter rmx = 10^ MeV (solid curves). For 
comparison, the corresponding ratios are also shown for the sudden decoupling ap- 
proximation [3] (dashed curves). 

FIG. 3. The change in the primordial "^Hc abundance, AY, as a function of decaying 
particle lifetime r for a decaying particle with a density parameter mix = 10^ MeV 
and three different baryon-to-photon ratios: f] = 10~^° (dashed curve), rj = 10~^'^ 
(solid curve) and rj = 10'^ (dotted curve). 

FIG. 4. The change in the primordial ^He abundance, Ay, for rj — 10~^-^, as a 
function of decaying particle lifetime r for a decaying particle with a density parameter 
rmx = 10^ MeV (solid curve) compared with the results of the sudden decoupling 
approximation [3] (dashed curve) and the sudden decoupling approximation with the 
interaction rates updated to correspond to the ones used here (dotted curve). 
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